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Abstract 



Single-molecule experiments in which force is apphed to DNA or RNA molecules have enabled 
important discoveries of nucleic acid properties and nucleic acid-enzyme interactions. These 
experiments rely on a model of the polymer force-extension behavior to calibrate the experi- 
ments; typically the experiments use the worm-like chain (WLC) theory for double-stranded 
DNA and RNA. This theory agrees well with experiments for long molecules. Recent single- 
molecule experiments have used shorter molecules, with contour lengths in the range of 1-10 
persistence lengths. Most WLC theory calculations to date have assumed infinite molecule 
lengths, and do not agree well with experiments on shorter chains. Key physical effects that 
become important when shorter molecules are used include (i) boundary conditions which 
constrain the allowed fluctuations at the ends of the molecule and (m) rotational fluctuations 
of the bead to which the polymer is attached, which change the apparent extension of the 
molecule. We describe the finite worm-hke chain (FWLC) theory, which takes into account 
these effects. We show the FWLC predictions diverge from the classic WLC solution for 
molecules with contour lengths a few times the persistence length. Thus the FWLC will al- 
low more accurate experimental calibration for relatively short molecules, facilitating future 
discoveries in single-molecule force microscopy. 

Key words: worm-like chain; single-molecule experiments; DNA; RNA; force-extension 
measurements; short chains; molecular stretching experiments; persistence length; contour 
length 
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Introduction 

Single-molecule force microscopy, in which force is applied to DNA or RNA molecules one at 
a time, has enabled important discoveries of the mechanical behavior of nucleic acid (NA) 
molecules and NA-enzyme interactions (P, 0). These experiments are usually calibrated 
and interpreted using the worm-like chain theory (WLC). The WLC predicts the average 
end-to-end extension 2; of a semiflexible polymer, given the force F applied to the ends of 
the chain and the values of two constant parameters. The first parameter is the molecule's 
unstressed total contour length L, which is proportional to the number of base pairs in a 
double-stranded nucleic acid polymer. The second parameter, the effective bending stiffness 
parameter or persistence length A, describes the molecule's local elastic behavior. Hence we 
expect the value of A to depend on the molecule type and the nature of the surrounding 
solution, but not on the contour length. (Additional parameters enter when the stretching 
force exceeds about 20 pN, or when the molecule is torsionally constrained.) The WLC has 
been successfully applied to long molecules of DNA (j^, 0) and RNA (j^ |^ and it is now the 
standard model used in calibrating and interpreting SM force microscopy (0). 

Existing treatments of the WLC theory do not agree well with experiments when the con- 
tour lengths of the NA molecules are short. In particular, experiments on dsDNA molecules 
~ 1-10 yield fit values of the persistence length which are physically unrealistic 
These apparent values differ from the accepted value, A ^ 50 nm, by factors of 
2-5. However, the persistence length is a material parameter, independent of the contour 
length. Therefore this apparent L-dependent persistence length amounts to a failure of the 
WLC, or at least its usual mathematical treatment, which assumes that L ^ A. In ad- 
dition, the classic WLC solution neglects some physical effects present in the experiments, 
including the bead(s) attached to the end(s) of the molecule and surface effects. In this 
paper we extend the WLC to include boundary conditions at the ends of the polymer and 
rotational fluctuations of the bead(s) attached to the molecule, and present solutions valid 
at finite L. Our "finite worm-like chain" (FWLC) solution gives more accurate predictions 
than the classic WLC solution for polymers with contour lengths less than a few times the 
persistence length. In related recent work, Kulic et al. have studied finite-length effects in 
the stretching of DNA with a kink( l. lf ) ; their results apply in the high-force limit. At the 
other extreme, references |l^ and ([isl i give numerical results for the "tethered particle", or 
zero-force, limit. 




Need for extensions of the WLC 

The interactions between proteins and nucleic acids are essential to cells. Processes such 
as transcription and translation, splicing, DNA copying, genome maintenance and DNA 
repair require NA-protein interactions. Premature aging, cancer, and even death can result 
from defects in NA-protein interactions. (We use the abbreviation NA to refer to single- or 
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double-stranded DNA, RNA or DNA-RNA hybrid molecules). Single-molecule experiments 
have given new insight into NA-protein interactions, yielding information that is difficult 
to determine by other methods (|l|). In single- molecule force microscopy, force is applied to 
assess NA mechanical behavior or NA-protein interactions ()l4( ). Typically one end of an NA 
molecule is chemically attached to a surface and the other end is attached to a bead (figure 
H)). Manipulation of the bead, usually by optical (flsh or magnetic f)l6h means, stretches the 
NA polymer. (In variants of this setup, the NA may be attached to beads at both ends, or an 
NA-binding protein may be attached to a surface.) Changes in the NA molecule extension 
and force can be measured or applied. 

A theory of NA stretching behavior is required to calibrate single-molecule force exper- 
iments. The average end-to-end extension 2; of a polymer depends on the applied force F. 
Using an accurate model of the force-extension behavior means that once the calibration is 
performed, knowledge of any two of the extension, force, and contour length determines the 
other. Thus, changes in experimentally measured extension at fixed force can be directly 
related to changes in contour length, for example when a processive enzyme "reels in" NA 
as it moves. 

In the absence of an accurate model of NA force-extension behavior, single-molecule force 
microscopy encounters several problems. First, experimentalists lack a check on whether the 
trap calibration is correct. Second, without an accurate force-extension theory, the contour 
length cannot be deduced from a measurement of the extension at fixed force. Therefore 
motion of motors which change the length of the polymer in time cannot be accurately 
measured. Third, force-extension measurements are frequently used to assess the number 
of polymers attached to a given bead. The attachment of polymers to the surface/bead 
occurs stochastically and can give more tethering molecules than desired. Measurement 
of the force-extension behavior is used to conffim the correct number of tether molecules. 
Fourth, experiments which study forced unfolding of secondary structure (jl, |^ [l3) require 
force-extension calibration to determine the free energy of the folded molecule. 

The worm-like chain theory (WLC) of NA elasticity has become the standard model for 
calibrating experiments. The WLC assumes that the polymer is an infinitely long, inexten- 
sible, isotropic rod with bending rigidity. The usual WLC solution gives the relationship 
between fra;tional extension z/L and force in the hmit L ^ 00 fl i However, phys- 
ical effects neglected in the usual treatment of the WLC are important for relatively short 
molecules, such as those used in several experiments (jl, |^ [l^, ITtI llih. 

Experiments have fit data to the usual solution to the WLC even for relatively short 
chains, because no alternative treatment has been available 1^0). The result is a fit value 
of the persistence length which depends on the contour length of the molecule, with fit values 
as low as 10 nm (Jq^). For double-stranded DNA, the accepted value of the persistence length 
is approximately 50 nm (depending on solution conditions ([lil lioi)). Experimental force- 
extension data, when fit to a theory that correctly incorporates finite-length effects, should 
recover the intrinsic value of A, independent of L. 
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Finite-length effects can be incorporated into tlie WLC tlieory by considering boundary 



conditions at tlie end of tlie cliain (|llL |2l|). Tlie constraints imposed by tlie boundary 



conditions alter the force-extension behavior. A second important physical effect typically 
neglected in the WLC theory is the rotational motion of the bead(s) attached to the molecule 
(figure EI). Experiments measure the position of the center of the bead, and typically one 
estimates the molecule extension by subtracting the bead radius. This estimate is correct 
only when the polymer-bead attachment is in line with the direction of the force {6 = 
in figure El) and the bead undergoes no rotational fluctuations. For long molecules (with 
contour lengths of tens of /xm), the fractional error introduced by subtracting the bead 
radius (typically 50-500 nm) can be neglected. However, for polymers with L of few hundred 
nanometers, the error caused by ignoring bead rotational fluctuations can be significant. 

Several other physical effects are present in experiments, but neglected in most treatments 
of the WLC. These include (i) the exclusion interaction between the molecule and the surface 
(which prevents the polymer from entering the solid surface), (m) the exclusion interaction 
between the bead and the surface, and {in) the exclusion interaction between the molecule 
and the bead. These interactions are most important for zero or very low applied force, and 
have been addressed in recent work (0, 0, |22 ) . Other interactions between the surface and 
the chain or bead may occur, such as electrostatic or van der Waals effects. In this paper, 
we will neglect these physical effects, which are important only for very low forces or very 
long polymers. Our goal in this paper is to present a unified model which includes the effects 
most important for single-molecule force microscopy experiments on dsNA molecules with 
~ 1-10; for dsDNA this means contour lengths of 50-500 nm. 



Theory 

The WLC model assumes an inextensible polymer: the contour length L of the molecule (the 
total length of the molecular backbone) cannot be changed by applied force. For double- 
stranded DNA, the overstretching transition which occurs for applied force around 60 pN 
is a dramatic change in the organization of the molecule; however, lower applied forces (up 



to approximately 20 pN) respect the inextensibility constraint (|23l . |2J). The polymer is 
assumed to possess an isotropic bending rigidity, characterized by the persistence length, 
A. The persistence length is the length scale over which thermal fluctuations randomize the 
chain orientation. (In principle, the rod also resists twist. However, in many experiments 
the twist is unconstrained and can be ignored.) We note that when considering single finite- 
length molecules, we are not working in any "thermodynamic limit". Therefore, different 
ensembles are „ot guaranteed to be equwatent 0HB. In this paper we work in the 
ensemble relevant to most experiments, where the applied force is fixed and calculate the 
extension. 

The chain energy includes the bending energy (Hooke's law in the chain curvature) and 
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the work done by the apphed force. The energy (in units of the thermal energy ksT) is 



E = -Fz + 



ds —K^ 







Here F (the force divided by the thermal energy ksT) is applied in the z direction, s 
denotes arc length, and the total extension of the chain is z. The curvature can be defined in 
terms of arc- length derivatives of the chain coordinate (figure |2|). If the chain conformation 
is described by a space curve r(s) and the unit vector tangent to the chain is t(s), then 

. Note that the chain extension can be calculated as 2; = jj" ds z -i. 









di 








ds 


We resca 


le I 



E 



ds' 



J2 



fz-t 



(2) 



where £ = L/A, f = FA, s' = s/A, and k' = kA. We will drop the primes in the remainder 
of the paper. 

To determine the extension for a given applied force requires averaging over different 
polymer conformations. This leads to a path integral formulation of the statistical-mechanics 
problem in the tangent vector to the chain (0, H^j). If the ends of the chain are held at fixed 
orientations, the partition function of the chain is 



Dt exp 



ds{-{d,ty-fz-t 



(3) 



where the integral in Dt is over all possible paths between the two endpoints of the chain 
with the specified orientations. The partition function can be rewritten as a propagator, 
which connects the probability distribution for the tangent vector at point s, ^(t,s) to the 
same probability distribution at point s': 



7/;(t, s) = j di' Z{i, s] t', s') tp{i', s' 



(4) 



From this relation, one can derive a Schrodinger-like equation which describes the s evolution 
o^^P i): 



ds 



+ / cos 9 ip. 



(5) 



Here is the two-dimensional Laplacian on the surface of the unit sphere and cos 6* = z • t. 
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Additional physical effects included in the FWLC 

Boundary conditions at ends of chain 

For relatively short NA molecules, the boundary conditions at the ends of the chain must 
be considered Elh . The boundary conditions are specified by two probability density 
functions, ip{t,s = 0) and ip{t,s = L). These enter the full partition function via 

Ztot = I dii dtf i,{ti, 0)Z(t„ 0; ty, L)ij{tf, L). (6) 



Bead rotational fluctuations 

Rotational fiuctuations of the bead can be explicitly included in the theory. Here we describe 
the extension due to a bead at one end of the chain, with the other end of the chain held at 
fixed orientation (figure^; the generalization to consider beads at both ends of the polymer 
is straightforward). In this case, the theory predicts the extension to the center of the bead, 
not just the molecule extension. We assume that the bead is spherical and the fiuctuations 
in the polymer-bead attachment are azimuthally symmetric about the z axis (figure El)- 

The unit vector fi points from the center of the bead to the polymer-bead attachment. 
Thus the bead's contribution to the energy of a fiuctuation with a given n is FR\\ ■ z, where 
R is the bead radius. Note that the energy is minimized when n = — z, and for other bead 
orientations the energy increases. This gives a restoring torque on the bead which increases 
as FR increases. We write the rescaled energy including this term as 



E = frh • z + / ds [ /z ■ t 



K 



(7) 



where r = R/A. The partition function of the chain is 



Dt exp 



- / ds 



-(d^if-fz-i 



constr 



(in exp (— /rn • z) 



where the first integral is over all possible paths between the two endpoints of the chain, and 
the second integral is over all allowable vectors fi. Below we evaluate the second integral, 
thereby performing the average over fiuctuations in n. 



Methods 

Here we describe the calculation of the force-extension relation, the main quantity of interest 
in single-molecule experiments. We must first calculate the tangent-vector probability dis- 
tribution ip{t, s) for all s along the chain. The distribution satisfies equation ^ above. This 
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PDE can be solved using separation of variables in s and t, where the angular dependence 
is expanded in spherical harmonics Q). 

oo 

^{i,s) = Y,^,{s)Y,o{i). (9) 

j=0 

We assume azimuthal symmetry, so only the m = spherical harmonics, with no depen- 
dence, are included. In the basis of spherical harmonics, the operator in equation (jSj) is a 
symmetric tridiagonal matrix H with diagonal terms 

H,., = (10) 

and off-diagonal terms 

= (2j + ^)(2j + 3)- ^''^ 
The vector of coefficients at point s is given by the matrix exponential of H: 

= e^^*(0). (12) 

This expression allows us to compute the probability distribution of the tangent vector 
orientation at any point along the chain. This result is exact if the infinite series of spherical 
harmonics is used. In practice, the series must be truncated for numerical calculations. Our 
calculations use = 30 unless otherwise specified. 

The partition function is then calculated from the inner product 

Z = = £)e^^*(s = 0), (13) 

= J2'^^{s = i)[e'%k'^kis = 0). (14) 

The extension at a given applied force is 

z IdlnZ 



L i df 

This formula applies for a chain of any length. 



(15) 



Finite-length correction 

To calculate the force-extension relation we must determine M = e^^ [equations ()12|1 and 
flT^ ]. Because H includes at least one positive eigenvalue, the entries M grows rapidly with 
i. This increase can lead to numerical overflow errors when computing M. However, we are 
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interested not in the entries of the matrix but in the logarithmic derivative of the partition 
function, and a rescahng can avoid the overflow problem. 

Let A = and denote by A* the largest eigenvalue of A. Then M = A^. If we define 
^ = y4/A*, we have M = XlA^. Thus A has eigenvalues with magnitude less than or equal 
to 1. The partition function can be written 

Z = Xi^^is = e)A'<if{s = 0). (16) 

The logarithm of the partition function is then 

In Z = £ln A, + ln[*^(s = = 0)]. (17) 

In the usual WLC solution, only the first term is considered |3), an approximation which 
is exact in the limit i ^ 1. The second term is the correction to InZ due to finite- length 
effects. 



Boundary conditions 

The boundary conditions at s = and s = i affect the force-extension relation, because 
they affect the partition function [equation (fT^ ] . The functional form of a specific boundary 
condition is given by the vector of coefficients ^'(s = 0), which is determined by the projec- 
tion of ijj{t, s) onto spherical harmonics. To apply different boundary conditio ns, we simply 
determine the partition function for different vectors ^{s = 0) and ^{s = i) (|2l[). 

In this paper, we consider three types of boundary conditions (figure |3I). First, in the 
"unconstrained" boundary condition the tangent vector at the end of the chain is free to point 
in any direction on the sphere (in Air of solid angle, figure |3^). In real experiments, however, 
the boundary conditions are more constrained. If the polymer is attached to a surface about 
a freely rotating attachment point, we might expect "half-constrained" boundary conditions 
(figure IDd), where the tangent vector at the end of the chain can point in any direction on 
the hemisphere outside the impenetrable surface. Many experiments appear to implement 
half-constrained boun dary conditions, because they use a flexible attachment between the 
chain and the surface (|l3h. We also consider perpendicular boundary conditions, where the 
tangent vector at the end of the chain is parallel to the z axis, normal to the surface (figure 

Et). 

For the unconstrained boundary condition, ip(t) is independent of cos 6'. Therefore ^ = 
(1, 0, ■ ■ ■ ,0). For the half-constrained boundary condition, the leading coefficients of ^ are 
(1, 0.8660, 0, -0.3307, 0, 0.2073, 0). Finally, for the perpendicular boundary condition the 
coefficients of ^ are all equal to 1. (Note that for the computation of the force-extension 
relation, it is not necessary to properly normalize the probability distribution, because we 
are computing the derivative of the logarithm of Z. Our expressions for the probability 
distribution vectors will neglect the constant normalization factor.) 
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Bead rotational fluctuations 

The bead rotational fluctuations lead to an effective boundary condition (at the end of the 
polymer) that depends on applied force and bead radius. For the single-bead experimental 
geometry (figure H^), equation (jHI) gives the partition function of the system including bead 
rotational fluctuations. We can perform the integral over fi to find the probability distri- 
bution of t(^), that is, is the effective boundary condition at the end of the polymer. The 
integral over n is 

g{m,i)= I rfne-^'-"^ (18) 
-/constr 

Because the direction of n is constrained relative to the chain tangent t(^), this integral is a 
function of t(£). By azimuthal symmetry, it depends only on the scalar t(£) ■ z. Accordingly 
we can express equation (fTHj) in the form g{t{t) -z), where g is the probability distribution of 
tangent angles at s = (Note that we neglect the normalization constant for g.) Expanding 
in spherical harmonics, we write 

oo 

gm ■ z) = ^?(cos7) = J2^A^ = i)Y,o{l). (19) 

j=0 

The effective boundary condition depends on both the applied force and the radius of the 
bead. The physical character of the chain-bead attachment determines the constraints in 
the integral of equation (fTH|l. and therefore controls the 

For the two-bead experimental geometry, the partition function includes integrals over 
bead rotational fluctuations for both ends of the chain, and the effective boundary condition 
applies at both ends of the polymer. 

Half-constrained boundary conditions 

Suppose that the polymer-bead attachment is half-constrained, so the tangent vector is free 
in a hemisphere (figure EJd). Then the integral of equation (|18|) is constrained by t(£) • n < 0, 
or 

^(t • z) = / dn e-f'^\ (20) 

Jtn<0 

To evalute the integral, we choose a polar coordinate system where t • z = cos 7, t points 
along the polar axis, and = corresponds to the z direction. Therefore z = (sin 7, 0, cos 7), 
n = (sin 9 cos 0, sin 9 sin 0, cos 9) , and n • z = sin 9 cos sin 7 + cos 9 cos 7. We can then write 
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g in terms of the Bessel function Jq. 



d<p j dcosO e--^^(^''^^™^'^'^''^^+=°'^^™'^^), (21) 



1 





2tt 



27T j dcosO e-^''^°''^'°'"'Jo{-ifrsmesmj). 



(23) 



Note that exp{zcosx) = TrJo{iz), and the integral from to 27i is even about vr. We show 
the values of this function in figure Elfor different apphed forces and bead radii. As expected, 
for small applied force and small bead radius, the probability distribution approaches a 
constant value (independent of cos 7). However, for large applied force, the probability 
distribution approaches the half-constrained distribution one would expect in the absence of 
bead rotational fluctuations. 



Numerical methods 

Because cos^ is integrated over negative values in equation (j^Hj) . the term g-/'^ cose cos 7 
verges as fr increases. However, we only calculate the dependence of g on cos 7 (correct 
normalization is not required). Therefore we replace the term g-/^ cos cos 7 ^j^g integral 
with the term e-/^'{™"'^™«T+i). 

To determine the effective boundary condition due to the rotationally fluctuating bead 
requires that we expand integrals of the form given in equation (PHj) in spherical harmonics. 
Direct numerical projection of the integral onto spherical harmonics leads to large errors, be- 
cause numerical integration of rapidly oscillating functions (such as the higher-order spherical 
harmonics) is inaccurate. To solve this problem, we used an interpolating basis (reference 
f|3o[ ). section 3.1.4). This allows the coefficients of the spherical harmonics to be determined 
by evaluation of equation at specific points (corresponding to the zeros of the Legendre 
polynomials). The integral in equation (|^ was evaluated numerically using Gauss- Legendre 
quadrature. 



Results 

Here we delineate the regimes where the usual WLC solution applies, and where the FWLC 
is required for good agreement with experiment. The classic WLC predictions becomes less 
accurate as the contour length of the chain decreases. The differences between the WLC and 
FWLC also depend on applied force, boundary conditions, and bead radius. 

As expected, the FWLC predictions converge to the WLC predictions as the polymer 
contour length increases. Figure IHl shows the predicted fractional extension z/L as a function 
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of the chain contour length L (for fixed apphed force) . Predictions calculated from the classic 
WLC solution are independent of contour length. However, for the FWLC the predicted 
fractional extension deviates from the infinite value for smaller values of L. In all cases, the 
FWLC predictions converge to the classic WLC results as L ^ oo. 

Error threshold contour length 

To summarize the different predictions of the FWLC and WLC, we calculate the "error 
threshold contour length" Lg. As shown in figureEl the FWLC and WLC predictions diverge 
as L is decreased. We decrease L and compare the FWLC and WLC predictions. The error 
threshold contour length is defined as the contour length where the relative difference 
between the two predictions first reaches 5%. For contour lengths L < L^, the FWLC model 
is necessary for accurate predictions. By characterizing the dependence of Lg on boundary 
conditions, applied force, and bead radius, we can understand which parameters control the 
difference between FWLC and WLC predictions. 

Figure [7| shows the error threshold length as a function of applied force for different 
boundary conditions and bead radii. For all conditions described, Le is largest for low applied 
force, where the predicted extension is most altered by finite-length effects. This behavior is 
intuitively reasonable. When the applied force tends to infinity, only one chain conformation 
is possible — the chain is perfectly straight, aligned in the direction of the force. In this case, 
modifying the theory to include boundary conditions or bead rotational fluctuations does 
not alter the polymer conformation. The characteristic propagation length of deformations 
in a classical elastic rod (that is, neglecting thermal fluctuations) is \/ AjF . As the force 
increases, the effects of any constraints on the ends of the polymer have decreasing effect 
on the conformation. Similarly, when the applied force is large bead rotational fluctuations 
require a large amount of energy and are suppressed. 

When the applied force is low, many different chain conformations are probable. There- 
fore the boundary conditions and bead rotational fluctuations exert greater influence on 
the chain extension. Unconstrained boundary conditions maximize the number of allowable 
conformations. In this case the entropy is maximized, and the polymer extension is low: 
increasing the extension requires excluding many chain conformations, which requires more 
work. More constrained boundary conditions restrict the number of conformations that are 
possible, leading to lower entropy and therefore a larger extension. At low force, the most 
constrained boundary conditions typically lead to the largest extension. We illustrate this 
idea with the limit of a perfectly rigid rod which is constrained at one end to be perpendic- 
ular to the surface. The infinite rigidity means that the molecule must be perfectly straight. 
The boundary conditions require that the polymer take only one conformation, perpendic- 
ular to the surface. In other words, the boundary conditions severely restrict the number 
of available conformations, and thereby increase the predicted extension. (If the boundary 
conditions are unconstrained, the rod rotates to point in many different directions and its 
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average extension is lower). In addition, for low applied force bead rotational fluctuations 
are larger, leading to a larger alteration of the extension due to bead fluctuations. 

Force-extension behavior 

We illustrate the force-extension curves for different boundary conditions and contour lengths: 
L/A = 4 in figure IHl and L/A = 10 in figure IHl As described above, (^) the fractional dif- 
ference between the infinite and finite predictions is largest for low contour lengths and low 
applied forces, (u) the more constrained boundary conditions lead to the largest changes 
relative to the infinite model, and [iii) the more constrained boundary conditions typically 
lead to larger extensions at low force. Therefore, the force-extension curves for the FWLC 
with unconstrained boundary conditions are closest to the classic WLC predictions (figure 

El). 

Including bead rotational fluctuations in the model leads to an effective boundary con- 
dition that depends on applied force and bead radius, as shown in figure El The variation of 
the effective boundary condition with force alters the shape of the force-extension curve. For 
high applied force, the fluctuations require more energy and so large-angle fluctuations oc- 
cur less frequently. Therefore, in the limit FR —>■ oo the probability distribution approaches 
that of the chain-bead attachment. However, for smaller FR, the bead rotational fluctua- 
tions "smear out" the probability distribution, making the effective tangent angle boundary 
condition approach a constant, independent of cos 7. For smaller FR, the the predicted 
extension is similar to that expected for no bead and unconstrained boundary conditions. 
This typically leads to a smaller extension than would occur in the absence of a bead. For 
larger FR, the effective boundary condition approaches the boundary condition that would 
occur with no bead. In this case, there is little difference between the predicted extension 
with and without a bead. 

We note that two opposing trends occur as the bead radius varies: when the bead is 
smaller, the fractional error made by subtracting R from the measured extension to estimate 
the molecule extension is smaller, simply because one is subtracting a smaller value. However, 
for smaller R the restoring torque on the bead for a given F is also smaller, which leads to 
larger angular fluctuations for smaller beads. 

At very low forces our FWLC theory gives inaccurate predictions, because it does not 
explicitly include the exclusion interaction between the bead and the wall. If the force is 
very low, the predicted extension would imply that the bead overlaps with the wall (or the 
other bead). Typically we find that applied force > 0.08 pN is required to prevent predicted 
bead-wall overlap in our theory. Although the FWLC does lead to unphysical predictions 
at very low forces, the simple model of the bead rotational fluctuations is useful for larger 
values of the force. 
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Fitting FWLC curves to the classic WLC solution 

Here we illustrate the typical errors that occur if one fits force-extension curves generated 
by the FWLC to the usual WLC solution. We Erenerated force-extension data with the 
FWLC, and fit to the WLC solution of Bouchiat et al (7). We find that the contour length 
is typically fit well by the WLC model, but large errors in the fitted value of the persistence 
length occur. The apparent persistence length decreases with contour length f figure ITUI) . as 
observed experimentally in the stretching of short DNA molecules (0). The magnitude of 
the error increases as the boundary conditions become more constrained, or when beads are 
included in the model. 

We emphasize that our fitting result is not directly comparable to the fitting of experi- 
mental data, because our simulated "data" are not noisy. In fits of real data, the quality of 
the fit is affected by the amount of noise. However, our result does show the types of errors 
that occur: fitting "data" generated by the FWLC to the usual WLC solution leads to large 
errors in the persistence length for small L. By contrast, fitting to the FWLC recovers the 
correct value of A . 



Discussion 

In this paper we have developed the finite worm-like chain theory (FWLC), which predicts 
the force-elasticity behavior of polymers with short contour lengths (£ = L/A ~ 1-10). In 
addition to retaining contributions from the subleading terms of equation (fT7|) . the FWLC 
includes two physical effects neglected in the usual WLC treatment: (z) chain-end boundary 
conditions and (m) bead rotational fluctuations. Together, these effects may explain the 
apparent decrease in polymer persistence length which has been experimentally measured 
for short chains (0). 

A key result of this paper is the delineation of the regimes where we expect the usual 
WLC solution to fail. We demonstrate that the FWLC converges to the WLC in the limit of 
long contour length for all forces, boundary conditions, and bead radii. However, for shorter 
molecules the contour length crossover where the usual WLC solution becomes inaccurate 
depends on applied force, chain-end boundary conditions, and bead radius. For the FWLC 
and WLC to agree within 5% for a force-extension measurement with applied force of 0.1-10 
pN, requires L/A> 100 (L > 5000 nm for dsDNA). 

Although the FWLC improves on the usual WLC solution for short contour lengths, 
it does not include all physical effects which occur in single-molecule force microscopy. In 
particular, effects omitted from the FWLC are important for very low or zero applied force 
(such as in tethered particle motion experiments). These effects include the chain- wall, chain- 
bead, and bead-wall exclusion interactions. All of these interactions are most important at 
zero applied force. In our calculations, the average position of the bead is "outside" the wall 
for forces above 0.1 pN. For accurate predictions at lower applied force, exclusion effects 
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need to be included. Comparison of the FWLC with models which include the exclusion 
interactions can further define the regimes of applicability of this model. In the future, an 
improved theory which also applies at very low force may become necessary. 

The FWLC will be useful for experiments that use short NA polymers in single-molecule 
force microscopy. A theory that is correct for small L/A helps ensure that the experimental 
force calibration is correct, that least-squares fitting indeed recovers the correct values of the 
contour and persistence lengths, and that the number of polymers attached between surface 
and bead can accurately be determined. The FWLC will facilitate experimental work with 
shorter polymers and contribute to future discoveries in nucleic acid-enzyme interactions. 
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Figure Legends 

Figure [H 

Sketch of typical experimental geometries for single-molecule force microscopy with nucleic 
acid polymers. The total molecule extension is z. (a) One-bead geometry. The molecule is 
attached at one end to a surface (glass slide or coverslip) and at other end to a bead. Force 
is applied to the center of the bead, (b) Two-bead geometry. The molecule is attached at 
both ends to beads. Force is effectively applied to the centers of both beads. 

Figure [21 

Variables used to describe chain conformations. Force is applied in the z direction. The arc 
length along the chain is s, and the vector r(s) = {x{s),y{s), z{s)) is the coordinate of the 
chain at arc length s. The vector t(s) is the unit vector tangent to the chain at s. 

Figure 1^1 

Variables used to describe bead rotational fluctuations. The vector n points from the center 
of the bead to the polymer-bead attachment. The angle 6 between the z direction and the 
bead-chain attachment is defined by cos^^ = — n ■ z. The bead has radius R. 

Figure HI 

Sketch of the tangent-angle boundary conditions, (a) Unconstrained boundary conditions. 
The tangent vector at the end of the chain (open arrowhead) can point in all directions 
relative to the z axis (notched arrowhead) with equal probability, (b) Half-constrained 
boundary conditions. Due to a constraint (such as the presence of a solid surface normal 
to the z axis) the tangent vector at the end of the chain can point only in the upper half- 
sphere, (c) Perpendicular boundary conditions. The tangent vector is constrained to point 
only normal to the surface, in the z direction. 

Figure [5l 

Effective chain-end boundary conditions induced by bead rotational fluctuations. The un- 
normalized probability density g{t-z) is shown versus cos 7 = t-z. Top, bead radius R = 100 
nm. Middle, bead radius R = 250 nm. Bottom, bead radius R = 500 nm. The different 
curves correspond to different values of the applied force. In this calculation the attachment 
of the chain to the bead is half-constrained: in the absence of bead rotational fluctuations, 
the probability density is zero for cos 7 < and constant for cos 7 > 0. This step function is 
approached for large force and large R. For smaller values of FR the probability density is 
smeared out, and approaches a constant value (independent of cos 7) for small FR. 
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Figure 1^1 

Convergence of the FWLC-predicted extension to tlie classic WLC prediction for large L. 
The predicted fractional extension z/L is shown as a function of the contour length L (at 
fixed applied force). The different curves correspond to different theoretical assumptions: 
(i) the classic WLC solution, calculated using the method of ref. (0), (m) the FWLC with 
unconstrained boundary conditions, {iii) the FWLC with half-constrained boundary con- 
ditions , and {iv) the FWLC with perpendicular boundary conditions (see figure 0}. The 
different panels show different applied forces: top, 0.01 pN; middle, 0.1 pN; bottom, 1 pN. 
The polymer persistence length is A = 50 nm. Note that the WLC prediction is independent 
of contour length, and the FWLC predictions converge to the WLC prediction as L increases. 
The convergence occurs more quickly for higher applied force (see text, note different y-axis 
scales for different panels). For low applied force, the more constrained boundary conditions 
lead to larger predicted fractional extension (see text). 

Figure [71 

Error threshold length as a function of force. The error threshold length is the contour 
length where the FWLC- and WLC-predicted extensions differ by 5%. For values of L < Lg, 
the classic WLC solution gives significant errors in the predicted extension. The polymer 
persistence length is A = 50 nm. The error threshold length is largest at low force. 

Top: no bead. The different curves correspond to unconstrained, half-constrained, and 
perpendicular boundary conditions. The more constrained boundary conditions have larger 
error threshold lengths at low force than the unconstrained boundary condition. Comparable 
predicted extensions by the classic WLC solution and the FWLC over the range of forces 
and boundary conditions shown requires L > nm {L/A > 200). 

Middle: one bead. The boundary conditions are half-constrained on both ends of the 
chain. We calculate the extension predicted to the center of the bead, then subtract the bead 
radius and compare to the extension predicted by the WLC. Larger bead radii lead to larger 
values of Lg at small forces, while for larger values of the applied force Lg is independent 
of the bead radius. Comparable predicted extensions by the classic WLC solution and the 
FWLC over the range of forces and boundary conditions shown requires L > 8 x 10^ nm 
{L/A > 1600). 

Bottom: two beads. The boundary conditions are half-constrained on both ends of the 
chain. We calculate the extension predicted between the centers of the beads, then subtract 
twice the bead radius and compare to the extension predicted by the WLC. Larger bead 
radii lead to larger values of at small forces, while for larger values of the applied force, 
Lg is independent of the bead radius. The values of are larger for the two-bead case 
than the one-bead case. Comparable predicted extensions by the classic WLC solution and 
the FWLC over the range of forces and boundary conditions shown requires L > 10^ nm 
{L/A > 2000). 
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Figure HI 

Force-extension curves predicted by the FWLC model for L = 200 nm {L/A = 4). 

Top: no bead. The different curves correspond to the classic WLC solution and the 
FWLC with unconstrained, half-constrained, and perpendicular boundary conditions. The 
predicted extension for all boundary conditions converges at high force. For low applied force, 
the more constrained boundary conditions lead to larger predicted extension (see text). 

Middle: one bead. The different curves correspond to bead radii of 100, 250, and 500 nm. 
The predicted molecule extension (after the bead radius has been subtracted) agrees with 
the no-bead prediction for large force, but decreases more rapidly as the force decreases. For 
low applied force (below 0.05 pN) the predicted extension becomes negative, because the 
bead-wall exclusion interaction is neglected in the model. In this regime the theory gives 
unphysical predictions. 

Bottom: two beads. The different curves correspond to bead radii of 100, 250, and 500 
nm. The predicted molecule extension (after twice the bead radius has been subtracted) 
agrees with the no-bead and one-bead predictions for large force. For low applied force 
(below 0.08 pN) the predicted extension becomes negative, because the bead-bead exclusion 
interaction is neglected in the model. In this regime the theory gives unphysical predictions. 

Figure ^ 

Force-extension curves predicted by the FWLC model for L = 500 nm {L/A = 10). Com- 
pared to the curves for L = 200 nm (figure |HI) the effects of boundary conditions and beads 
are less important. 

Top: no bead. The different curves correspond to the classic WLC solution and the 
FWLC with unconstrained, half-constrained, and perpendicular boundary conditions. 

Middle: one bead. The different curves correspond to bead radii of 100, 250, and 500 nm. 
For low applied force (below 0.02 pN) the predicted extension becomes negative, because 
the bead-wall exclusion interaction is neglected in the model. In this regime the theory gives 
unphysical predictions. 

Bottom: two beads. The different curves correspond to bead radii of 100, 250, and 
500 nm. For low applied force (below 0.06 pN) the predicted extension becomes negative, 
because the bead-bead exclusion interaction is neglected in the model. In this regime the 
theory gives unphysical predictions. 

Figure ITOl 

Results of fitting the FWLC model predictions to the classic WLC solution. The fit value of 
the persistence length is shown as a function of contour length L. The polymer persistence 
length is A = 50 nm. 
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Top: no bead. The different curves correspond to tfie unconstrained, fialf-constrained, 
and perpendicular boundary conditions. Tfie apparent A decreases most for tlie lialf- 
constrained boundary conditions. Tlie apparent persistence lengtli decreases by more tlian 
a factor of two as L decreases from 10"^ to 100 nm. 

Middle: one bead. The different curves correspond to bead radii of 100, 250, and 500 nm. 
The decrease in the apparent persistence length is similar to that shown in the top panel, 
half-constrained boundary conditions. 

Bottom: two beads. The different curves correspond to bead radii of 100, 250, and 500 
nm. The apparent persistence length decreases more quickly with contour length than for 
the one-bead case. 
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